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ABSTRACT 

For  about  two  decades,  Structural  Dynamics  and  Coupled  Systems  Department  of  ONERA  has  developed, 
validated  and  applied  different  methods  for  modeling  and  predicting  internal  and  external  noise  of 
complex  structures  in  the  low  (LF),  medium  (MF)  and  high  (HF)  frequency  ranges.  These  activities  were 
initialized  for  sound  radiated  from  submarine  structures  and  then  applied  to  aeronautic  and  space 
structures  as  helicopters,  planes  and  launchers,  for  internal  noise  prediction. 

For  the  LF  and  MF  domains,  these  methods  are  based  on  Finite  Element  modeling  for  both  structures  and 
internal  fluids  and  on  the  integral  method  for  unbounded  external  fluids  while  for  the  HF  domain,  they  are 
based  on  Statistical  Energy  Analysis  (SEA). 

•  For  the  LF  range,  the  classical  modal  approach  is  used  for  both  structures  and  internal  fluids. 
The  coupling  is  then  performed  on  the  “fluid-structure  ”  common  interface. 

•  For  the  MF  range,  a  specific  method  called  “Onera-MF  method”  has  been  developed.  This 
method  uses  a  direct  approach  for  the  coupling  between  internal  fluids  and  structures  and  it  was 
implemented  initially  in  an  “in-house  ”  version  of  ADINA  code.  Now  it  is  implemented  around 
Nastran  code. 

•  For  the  HF  range,  the  usual  SEA  method  has  been  reformulated  in  order  to  take  into  account  the 
coupling  of  external  and  internal  heavy  fluids  with  a  structure.  A  ONERA  code  has  been 
developed. 

A  lot  of  validations  and  applications  of  the  methods  and  associated  codes  were  performed  by  numerous 
comparisons  between  theory  and  measurements. 

•  The  first  validations  were  performed  by  numerical  methods  on  submarine  structures  and  on 
helicopters  in  LF  and  MF  domains. 

•  Now  these  methods  (mainly  SEA  approach)  are  used  to  predict  internal  noise  inside  -helicopters 
(from  mainly  Eurocopter)  and  -Airbus  passengers-cabins/cockpits  in  MF  and  HF  ranges  ;  for 
different  sources  of  excitation:  Turbulent  Boundary  Layer  extended  over  fuselage,  Air 
conditioning  inside  cockpits  for  planes  and  Noise  coming  from  Gearbox,  coming  from  main  rotor 
and  coming  from  engines  air  intakes  for  helicopters. 

A  review  of  the  previous  methods  and  their  main  validations  and  applications  for  aircraft  internal  noise 
are  presented  herein. 


Paper  presented  at  the  RTO  AVT  Symposium  on  Habitability  of  Combat  and  Transport  Vehicles:  Noise ,  Vibration  and 
Motion”,  held  in  Prague,  Czech  Republic,  4-7  October  2004,  and  published  in  RTO-MP-AVT-1 10. 
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ORGANIZATION 


1.0  INTRODUCTION 

In  this  paper,  we  will  present  first  a  summary  of  the  different  methods  used  at  ONERA  within  the 
framework  of  the  prediction  of  internal  and  external  noise  of  complex  vibrating  structures  in  the  three 
different  domains  of  frequency:  LF,  MF  and  HF  domains.  For  FF  range,  classical  modal  approach  is  used 
for  both  structures  and  internal  fluids.  For  MF  range,  a  specific  method  called  “  Oner a-MF -method”  has 
been  developed.  This  method  is  implemented  in  an  "in-house"  version  of  ADINA  code.  For  HF  range, 
usual  SEA  method  has  been  reformulated  in  order  to  take  into  account  strong  coupling  of  heavy 
unbounded  fluids  with  a  structure.  A  ONERA  code  has  also  been  developed.  Main  applications  of  these 
methods  to  internal  noise  prediction  of  aeronautical  structures  will  be  presented  finally. 


2.0  REVIEW  OF  METHODS  USED  AT  ONERA  FOR  VIBRO ACOUSTIC 
PROBLEMS 

2.1  Numerical  method  for  LF  domain 

FF  domain  can  be  defined  as  the  frequency  band  where  spectral  responses  of  coupled  system  {external 
fluid} -{structure} -{internal  fluids}  reveal  isolated  resonances,  indicating  the  presence  of  an  underlaying 
modal  structure  of  an  associated  conservative  system.  These  resonances  come  from  the  first  elastoacoustic 
eigenmodes  of  coupled  system  {structure} -{internal  fluids}.  This  method  is  fully  detailled  in  [1,  2]  and  it 
is  based  on  the  classical  modal  approach  which  builts  up  the  coupling  of  the  two  modal  bases  of  first 
eigenmodes  of  both  structures  and  internal  fluids  over  the  “fluid-structure”  common  interface.  It  is  used  at 
ONERA  mainly  for  external  vibroacoustic  problems  of  submarine  structures.  For  internal  vibroacoustic 
problems,  we  prefer  to  use  a  MF  method  which  can  be  extended  in  that  case  to  FF  domain. 

2.2  Numerical  method  for  MF  domain 

In  MF  domain,  the  main  method  used  at  ONERA  is  a  method  called  “Onera-MF  method ”  [3,  4].  This 
method  has  been  developed  specifically  at  DDSS  Department  and  it  can  be  applied  to  vibroacoustic 
systems  defined  by:  -  a  continuous  and  complex  structure  which  is  constituted  of  elastic  or  viscoelastic 
materials,  dissipative  and  unhomogeneous  and  which  can  be  usually  discretized  by  classical  Finite 
Element  method  (FEM),  -  one  or  several  viscous  or  unviscous  internal  fluids,  compressible,  gas  or  liquids, 
occupying  acoustic  cavities,  -  a  compressible  external  unviscous  fluid,  gas  or  liquid,  occupying  an 
unbounded  domain. 


2.2.1  Formulations  used  for  vibroacoustic  system 

Coupling  between  external  fluid  and  structure  is  treated  by  method  of  integral  equations  which  solves  a 
Neumann  problem  linked  to  Helmhotz  equation.  Internal  fluid  is  treated  by  solving  a  Neumann  internal 
problem  which  is  formulated  in  terms  of  velocity  potential  field  for  fluid.  Internal  fluid  is  discretized  by 
Finite  Element  method.  Finear,  viscoelastic  and  solid  structure  is  also  treated  by  Finite  Element  method. 

2.2.2  Characteristics  of  Medium-Frequency  domain 

MF  domain  is  the  frequency-domain  where  modal  densities  of  systems  are  not  constant.  They  can  be 
either  high  or  either  poor  when  frequency  is  varying  within  the  MF-band.  In  this  domain,  characteristics  of 
materials  can  vary  with  frequency  and  it  is  not  possible  to  find  a  unique  conservative  system  associated  to 
the  analyzed  system  for  a  large  range  of  frequencies.  Moreover,  Finite  Element  models  used  for  the 
resolution  in  this  frequency-band  must  be  built  up  using  a  large  number  of  nodes  and  consequently  a  large 
number  of  degrees-of-freedom  (DOF) ;  which  leads  to  handle  large  matrices  systems. 
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That  is  why  standard  methods  of  modal  reduction  usually  used  in  low-frequency  domain  are  not  accurate 
when  frequency  increases.  The  use  of  such  methods  requires  to  extract  a  very  high  number  of  eigenmodes 
of  the  vibroacoustic  system,  at  prohibitive  cost  and  without  assuming  convergence  of  iterative  algorithms. 
A  direct  frequency-by-frequency  method  is  not  either  envisaged  at  a  reasonable  numerical  cost  because  we 
have  to  inverse  for  each  discrete  frequency  a  very  large  linear  complex  matrices  system  having  a  very  high 
number  of  physical  D.O.F.s.  That  is  why  a  specific  method  adapted  to  MF  domain  has  been  developed. 

2.2.3  Onera-MF  method 

The  method  is  based  on  the  use  of  a  time  integration  and  Fourier  transform  technique  which  is 
summarized  below.  The  MF  broad  band  is  divided  in  several  narrower  frequency  sub-bands.  The  linear 
complex  matrices  system  is  then  solved  for  each  MF  sub-band  in  three  steps: 

1)  The  short  time  scale  which  is  associated  to  the  central  frequency  of  sub-band  is  analytically  treated  in 
frequency  domain  using  a  frequency  shift  technique. 

2)  The  long  time  scale  which  is  associated  to  the  bandwidth  is  solved  by  an  usual  direct  step-by-step 
numerical  integration  method  in  time  domain  (Newmark  method). 

3)  The  MF  solution  is  reconstructed  in  frequency  domain  by  Shannon  theorem  using  only  the  time 
solution  (at  sampling  time  steps)  constructed  in  step  2)  before. 

MF  method  has  led  to  develop  codes  which  allow  us  to  calculate  the  vibroacoustic  response  of  structures 
coupled  or  not  to  internal  or  external  fluids.  These  codes  were  initially  articulated  around  an  nin-housen 
version  of  ADINA  code.  Now,  MF  method  is  implemented  around  MSC/Nastran  code  within  I-DEAS 
environment. 

2.3  SEA  method  for  HF  domain 

Statistical  Energy  Analysis  (SEA)  is  adapted  to  HF  domain  which  can  be  defined  as  the  frequency  band 
where  spectral  responses  of  coupled  system  {external  fluid} -{structure} -{internal  fluids}  have  no 
resonances  nor  local  accidents.  In  this  domain,  dynamic  behaviors  are  very  smoothed,  indicating  that 
modal  density  of  conservative  system  associated  to  coupled  problem  {structure} -{internal  fluids}  is 
uniforme  in  frequency  and  is  high.  SEA  method  used  at  ONERA  and  its  associated  software  are  directly 
coming  from  a  reformulation  of  basic  SEA  method  [6].  This  reformulation  [7,  8]  was  performed  in  order: 

•  to  avoid  some  restrictive  hypotheses  of  the  basic  formulation  of  SEA  (iso-energy  of  resonant 
modes  of  a  subsystem,  iso-contribution  in  the  coupling  between  of  a  resonant  mode  of  a 
subsystem  to  any  resonant  mode  of  another  subsystem), 

•  to  improve  determination  of  new  coupling  loss  factors, 

•  to  introduce  a  model  for  estimating  the  envelope  of  SEA  prediction  around  averaged  response, 
due  to  uncertainties  of  SEA  parameters. 
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3.0  VALIDATIONS  AND  APPLICATIONS 

Validations  on  LF  method  have  been  performed  on  coupled  {structure} -{external  heavy  fluid}  [5,  1 1]. 

First  validations  of  “Onera-MF  method ”  were  focused  on  the  "own-noise"  of  sonar-domes  excited  by 
turbulence  boundary  layer  [10].  More  recently,  other  validations  concern  more  academic  systems:  a  cavity 
filled  with  water  and  coupled  to  a  plate  [18],  a  plate  coupled  to  external  air  and  excited  by  a  mechanical 
force  [12],  a  plate  coupled  to  external  air  and  excited  by  the  Turbulent  Boundary  Layer  [19],  a  cylinder  in 
air  and  in  water  [13].  These  three  latest  systems  also  allowed  to  validate  SEA  method  for  light  and  strong 
coupling  between  a  structure  and  unbounded  external  fluid.  SEA  method  has  also  been  validated  for  a  box 
which  is  representative  of  the  structure  of  an  helicopter  cabin  [16]. 

3.1  Applications  for  helicopters 

Main  application  of  MF  method  concerns  Dauphin  helicopter  of  Eurocopter  which  is  excited  by  random 
forces  due  to  Main  Gearbox  (MGB)  emergent  frequencies  and  applied  to  Transmission  deck  (see  Fig.  1) 
by  BTP  struts  [9].  Simplified  FE  model  used  for  MF  analysis  is  presented  at  Fig.  2  and  comparison 
between  in-flight  measurements  and  FE  computation  on  cabin  pressure  levels  obtained  for  main  MGB 
emergent  frequencies  is  presented  at  Fig.  3. 


I — I  ALU 

i — i  Honeycomb 
1 - 1  Sandwich 


Figure  1:  View  of  Dauphin  SA  365N  helicopter  and  selected  parts  for  FEM  analysis 
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Figure  2:  FE  structural  meshing  of  Dauphin 
fuselage  used  to  calculate  MF  responses 


Microphone  located 
inside  cabin-cavity 


Figure  3:  Comparison  between  numerical 
prediction  and  in-flight  measurement  of 
internal  pressure  levels  for  different  MGB 
emergent  frequencies 


Application  of  SEA  method  has  been  performed  on  internal  noise  prediction  within  Ecureuil  helicopter 
cabin  [15]  which  is  submitted  in-flight  to  four  main  different  sources  of  excitation  (presented  at  Fig.  4). 

Fig.  5  shows  comparison  between  the  envelope  of  SEA  prediction  and  in-flight  measurement  within  cabin 
for  this  helicopter. 
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Figure  4:  View  of  Ecureuil  helicopter 
subjected  to  different  sources  of  excitation 


Figure  5:  Comparison  between  SEA 
prediction  and  in-flight  measurement  of 
Ecureuil  cabin-cavity  induced  noise 
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ORGANIZATION 


3.2  Applications  for  planes 

Onera-MF  method  and  SEA  method  have  been  both  applied  to  the  prediction  of  internal  noise  inside 
cabin-passengers  of  a  Falcon  (from  Dassault)  airframe  [14,  17]  when  structure  is  excited  by  a  set  of  four 
mechanical  excitations  applied  at  left  engine  anchorage  points,  as  shown  in  Fig.  8.  Structure  of  fuselage 
and  cabin-passengers  acoustic  cavity  have  been  discretized  by  FE  method  (see  Fig.  9)  for  MF  approach 
and  all  the  structure  was  cut  up  in  nine  subsystems  (see  Fig.  10)  for  SEA  analysis. 


Figure  6:  Falcon  2000  test  airframe 


Figure  7:  Interior,  looking  to  the  back,  with 
microphones  in  place 


Figure  8:  Applied  mechanical 
excitations 


Figure  9:  Semi-longitudinal  section  of  Falcon  fuselage 
FE  meshing 
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Pressurization  partition 
Isolation  partition 


section  1 


section  2 

Length  tl 

Cavity  I 


section  3 


section  4 
Length  T5 

Cavity  2 


Mechanical  subsystems 


Acoustic  subsystems 


Lengths :  tl,T5 
Partitions  :  C4,  CP26,  C39 
Frames  :  C30.  C33 


Cavity  1 
Cavity  2 


Figure  10:  Breakdown  of  Falcon  airframe  into  SEA  subsystems 


Fig.  11  and  Fig.  12  illustrate  comparison  between  measurements  and  the  two  predictive  methods  on  the 
mean  quadratic  pressure  obtained  inside  cabin-passengers  cavity,  for  LF  and  MF  band  [0,  400Hz]  and  for 
HF  domain  [125,  1250Hz]. 


o  Front Y  excitation 

'fD  (N - 

Q_  1 


Frequency  (Hz) 


Figure  11:  Average  quadratic  pressure  for 
passengers-acoustic  cavity 


Frequency  (Hz) 


Figure  12:  Mean  quadratic  pressure  inside  the 
two  acoustic  cavities 


Another  main  application  of  SEA  method  concerns  prediction  of  TBL  excitation  induced  noise  inside 
cockpit  of  an  Airbus  A340  in-flight  (see  Fig.  13).  This  application  allowed  to  validate  a  SEA  model 
developed  at  ONERA  to  quantify  power  input  injected  to  structure  by  TBL  external  excitation  for  large 
structures. 
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Figure  13:  Comparison  between  SEA 
prediction  and  in-flight  measurement  of 
internal  noise  inside  A340-cockpit  cavity 
for  external  TBL  excitation 


4.0  CONCLUSIONS 

Some  efficient  methods  developed  at  ONERA  allowing  to  deal  with  external  and  internal  vibroacoustic 
problems  of  complex  structures  in  LF,  MF  and  HF  domains  and  their  applications  to  aeronautical 
structures  have  been  presented  in  this  paper.  Other  future  applications  of  SEA  method  are  planned  in  next 
few  months  for  two  large  aircraft:  helicopter  Super-Puma  of  Eurocopter  and  Airbus  A3  80. 
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Detailed  Analysis  or  Short  Description  of  the  AVT-110  contributions  and  Question/Reply 

The  Questions/ Answers  listed  in  the  next  paragraphs  (table)  are  limited  to  the  written  discussion  forms 
received  by  the  Technical  Evaluator.  The  answers  were  normally  given  by  the  first  mentioned  author- 

speaker. 


P23  J-M.  David  ‘Internal  Noise  Prediction  Within  Aircraft  Cabin’  (ONERA,  FR) 

A  third  contribution  to  this  very  important  problem,  taking  into  account  with  the  habitability  and  the 
resulting  safety  problems:  the  author  extensively  described  the  LF,  MF  (uneasy  to  define)  and  HF  (High 
Frequency)  approaches:  the  classical  MEF  under  100  Hz,  a  specific  ONERA-MF  method  (tested  a.o.  for  a 
cavity  filled  with  water)  implemented  in  the  Nastran  Code  ,  and  a  statistical  energy  analysis  (SEA)  over 
300  Hz.  Example  are  given  (Dauphin  helicopter  (MF),  Ecureuil  helicopter,  Falcon  (SEA),  ...The  author 
concluded  on  the  future  applications  of  SEA  method  planned  in  next  few  months  for  two  large  aircrafts: 
the  helicopter  Super-Puma  of  Eurocopter  and  the  Airbus  A3  80 

Discussor’s  name:  B.  Masure 

Q.  For  cavities  filled  with  water,  is  the  water  considered  as  a  compressible  medium? 

R.  It  doesn’t  matter  in  MF  method  that  the  water  is  considered  as  a  compressible  medium.  In  MF 
approach,  the  dynamics  of  both  structures  and  internal  fluids  is  included  within  the  formulation  and  also 
within  the  model.  Only,  in  modal  approach,  we  need  to  consider  the  water  as  an  incompressible  fluid.  In 
MF  method,  fluids  can  be  compressible  or  uncompressible. 

Q.  This  question  deals  with  the  excitation  forces  in  your  calculation  methods.  For  instance,  in  your 
presentation,  you  spoke  of  the  excitations  of  a  structure  associated  with  a  turbulent  boundary  layer  (TBL). 
Could  you  explain  how  exactly  you  do  to  define  these  excitations? 

R.  The  taking  into  account  of  excitation  by  TBL  depends  on  the  approach  used.  In  SEA  method,  a  model 
of  mean  power  input  due  to  TBL  excitation  was  developed  and  this  model  uses  the  mean  wall  pressure  and 
the  correlation  lengths.  In  MF  approach,  the  model  is  defined  by  the  data  of  a  cross-spectral  density 
function.  A  finite  element  discretisation  of  this  model  is  performed,  giving  us  a  spectral  matrix  of 
excitation. 

Discussor’s  name:  D.  Chan 

Q.  How  do  you  determine  the  boundary  between  the  MF  and  HF  frequency  domains? 

R.  We  don’t  have  any  criteria  to  determine  this  boundary.  We  are  able  to  determine  this  boundary  when 
we  have  measurements.  In  fact,  this  boundary  doesn’t  exist  really  but  that  we  are  sure  is  the  difference 
between  LF  modal  domain  and  HF  domain  and  this  intermediate  domain  is  the  MF  domain.  Only 
experience  can  define  the  MF  domain.  We  consider  HF  domain  as  the  domain  where  modal  densities  are 
high  and  constant  with  frequency.  If  it  is  not  the  case  for  a  complex  structure,  we  can  conclude  that  we 
still  are  in  the  MF  domain  and  the  boundary  between  MF  and  HF  is  not  yet  reached! 
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